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We report the observation of the b —* d penguin-dominated decay B° — ► K*°K*° with a sam- 
ple of 383.2 ± 4.2 million BB pairs collected with the BABAR detector at the PEP-II asymmetric- 
energy e + e _ collider at the Stanford Linear Accelerator Center. The measured branching fraction 
is B(B° -> K*°K*°) = [1.28t°;3o ±0.11] x 10~ 6 and the fraction of longitudinal polarization is 
f L (B° -> K*°K*°) = 0.80±o.i2 ±0.06. The first error quoted is statistical and the second sys- 
tematic. We also obtain an upper limit at the 90% confidence level on the branching fraction for 
J3(B° -> K*°K*°) < 0.41 x 10~ 6 . 
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PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh 



The study of the branching fractions and angular dis- 
tributions of B meson decays to hadronic final states 
without a charm quark probes the dynamics of both weak 
and strong interactions, and plays an important role in 
understanding CP violation. Decays proceeding via elec- 
troweak and gluonic b — > d penguin diagrams have only 
recently been measured in the decays B — > fry [l[ and 
B° -> K°K° 0. On the other hand, the charmless de- 
cay B° — > K *°K*° proceeds through both electroweak 
and gluonic b — > d penguin loops to two vector parti- 
cles (VV). The Standard Model (SM) suppressed decay 
B° — > K*°K*° could appear via an intermediate heavy 
boson. 

Theoretical models in the framework of QCD factor- 
ization predict the angular distribution of the VV de- 
cays of the B meson, as measured by the longitudinal 
polarization fraction Jl, to be ~ 0.9 for both tree- and 
penguin-dominated decays Q. However, recent measure- 
ments of the pure penguin VV decay B — > <f>K* indi- 
cate 0.5 Several attempts to understand this 
unexpected value of /l within or beyond the Standard 
Model have been made [5j. Further information about 
decays related by SU (3) symmetry may provide insight 
into this polarization puzzle and test factorization mod- 
els. A time-dependent angular analysis of B° — > K*°K*° 
can distinguish between penguin annihilation and rescat- 
tering as mechanisms for the value of /l observed in 
B -> cj)K* %. The B° -> K*°K*° mode can also be 
used within the SM framework to help constrain the an- 
gles a and 7 of the Unitarity Triangle 0- 

Theoretical calculations for B° — > K*°K*° branching 
fractions cover the range (0.16-0.96) x 10 -6 [§]. Recently, 
Beneke, Rohrer, and Yang @ predicted (O.eloiitoi) x 
10~ 6 and f L = 0.69 ± 0.01±g;^. Experimentally, upper 
limits on the branching fractions at the 90% confidence 
levelJC.L.) of 22 x 10~ 6 and 37 x 10~ 6 exist for B° -> 
K*°K*° and B Q -> K* Q K* Q , respectively 0. 

We report measurements of the branching fraction 
and the fraction of longitudinal polarization for the de- 
cay mode B° — > K*°K*°, with explicit consideration 
of non-resonant backgrounds and interference from K*° 
Tt'* (1430) . We place an upper limit on the branching 
fraction of B° — ► K*°K*°, where we use the notation 
K*°K*° to also represent K*°K*°. Charge-conjugate 
modes are implied throughout and we assume equal pro- 
duction rates of B + B~ and B a B° . 

This analysis is based on a data sample of 383.2 ± 4.2 
million BB pairs, corresponding to an integrated lumi- 
nosity of 348 fb , collected with the BABAR detector at 
the PEP-II asymmetric-energy e + e~ collider operated 
at the Stanford Linear Accelerator Center. The e + e~ 
center-of-mass (cm.) energy is ^fs = 10.58 GeV, cor- 
responding to the T(AS) resonance mass (on-resonance 



data). In addition, 36.6 fb~ of data collected 40 MeV 
below the T(45) resonance (off-resonance data) are used 
for background studies. 

The BaBAR detector is described in detail in Rcf. 



III. 



Charged particles are reconstructed as tracks with a 5- 
layer silicon vertex detector and a 40-layer drift chamber 
inside a 1.5-T solenoidal magnet. An electromagnetic 
calorimeter is used to identify electrons and photons. 
A ring-imaging Cherenkov detector is used to identify 
charged hadrons and provides additional electron identi- 
fication information. Muons are identified by an instru- 
mented magnetic-flux return. 

The B° -> K*°K*° and B° -> K*°K*° candidates 
are reconstructed through the decays K*° — > K + ir~ and 
K*° — > K~tt + . The differential decay rate, after inte- 
grating over the angle between the decay planes of the 
vector mesons, for which the acceptance is uniform, is 



d 2 T 



1 



r d cos 9id cos 02 



1-/. 



L ■ 2 

sin 



6>i sin 2 6> 2 +/ L cos 2 9t cos 2 6» 2 , 
(1) 

where 0\ and #2 are the helicity angles of the K* a or 
K*°. The helicity angle of the K*° (K*°) is defined as the 
angle between the K + (K~) momentum and the direction 
opposite to the B meson in the K*° (K*°) rest frame 12 1. 

The charged tracks from the K*° decays are required 
to have at least 12 hits in the drift chamber and a trans- 
verse momentum greater than 0.1 GeV/c. The tracks are 
identified as either pions or kaons by measurement of the 
energy loss in the tracking devices, the number of photons 
measured by the Cherenkov detector and the correspond- 
ing Cherenkov angles. These measurements are com- 
bined with calorimeter information to reject electrons, 
muons, and protons. We require the invariant mass of 
the K*° candidates to be 0.792 < m Kn < 1.025 GeV/c 2 . 
A B meson candidate is formed from two K*° candi- 
dates, with the constraint that the two K*° candidates 
originate from the interaction region. 

B meson candidates are characterized kinematically by 
the energy difference AE = Eg — \fsj2 and the energy- 

1/2 

substituted mass ibes = [( s /2 + Pi • ps) 2 /i? 2 — p%] , 
where (Ei, pi) and (Eb, Pb) are the four-momenta of the 
T(4S) and B meson candidate, respectively, and the as- 
terisk denotes the T(4S) rest frame. The total event 
sample is taken from the region —0.08 < AE < 0.2 GeV 
and 5.25 < tties < 5.29GeV/c 2 . Events outside the re- 
gion \AE\ < 0.07 GeV and 5.27 < m ES < 5.29 GeV/c 2 are 
used to characterize the background. The average num- 
ber of signal B meson candidates per selected data event 
is 1.03. A single B meson candidate per event is chosen 
as the one whose fitted decay vertex has the smallest x 2 - 
MC simulations show that up to 4% (1.6%) of longitudi- 
nally (transversely) polarized signal events are misrecon- 
structed, with one or more tracks originating from the 
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other B meson in the event. 

To reject the dominant background consisting of light- 
quark qq (q = u,d,s,c) continuum events, we require 
| cos 6t I < 0.8, where 9t is the angle, in the cm. frame, 
between the thrust axes [3] of the B meson and that 
formed from the other tracks and neutral clusters in the 
event. We create a Fisher discriminant T to be used 
in the maximum-likelihood (ML) fit, constructed from 
a linear combination of five variables: the polar angles 
of the B meson momentum vector and the B meson 
thrust axis with respect to the beam axis, the ratio of 
the second- and zeroth-order momentum-weighted Leg- 
endre polynomial moments of the energ y flo w around the 
B meson thrust axis in the cm. frame [14j . the flavor of 
the other B meson as reported by a multivariate tagging 
algorithm [HI , and the boost-corrected proper-time dif- 
ference between the decays of the two B mesons divided 
by its variance. The second B meson is formed by creat- 
ing a vertex from the remaining tracks that are consistent 
with originating from the interaction region. 

We suppress background from decays to charmed 
states by removing candidates that have decays consis- 
tent with D~ -> A+ 7r 7r and an invariant mass in the 
range 1.845 < m K +^- v - < 1.895 GeV/c 2 . We reduce 
backgrounds from B a — > cj)K* by assigning the kaon 
mass to the pion candidate and rejecting the event if 
the combined invariant mass of the two charged tracks 
is between 1.00 and 1.04GeV/c 2 . Finally, we require the 
cosine of the helicity angle of both A* candidates to 
be less than 0.98 to reduce continuum background and 
avoid the region where the reconstruction efficiency falls 
off rapidly. 

We use an extended unbinned ML fit to extract the sig- 
nal yield and polarization simultaneously for each mode. 
The extended likelihood function is 



(2) 



We define the likelihood £j for each event candidate i 
as the sum of rijVj{xi\ dj) over four hypotheses j (sig- 
nal, qq background, if* (1430) and BB backgrounds as 
discussed below), where Vj(xi \ dj) is the product of the 
probability density functions (PDFs) for hypothesis j 
evaluated for the i-th event's measured variables x il rij 
is the yield for hypothesis j, and N is the total num- 
ber of events in the sample. The quantities dj represent 
parameters in the expected distributions of the measured 
variables for each hypothesis j. Each discriminating vari- 
able Xi in the likelihood function is modeled with a PDF, 
where the parameters dj are extracted from MC simula- 
tion, off-resonance data, or (ibes, AE) sideband data. 

The seven variables Xi used in the fit are TOeSj AE, 
T , and the invariant masses and cosines of the helicity 
angle of the two A* candidates. Since the correlations 



among the fitted input variables are found to be on aver- 
age ~ 1% with a maximum of 4%, we take each Vj to be 
the product of the PDFs for the separate variables. The 
effect of neglecting correlations is evaluated by fitting 
ensembles of simulated experiments in which we embed 
signal and background events randomly extracted from 
fully-simulated MC samples. 

The two invariant mass and helicity angle distribu- 
tions for each A* meson are indistinguishable and so 
we use the same PDF parameters for both A* candi- 
dates. Peaking PDF distributions are described with an 
asymmetric Gaussian or a sum of two Gaussians. The 
transverse (longitudinal) helicity angle distributions are 
described with a cos 2 9 (sin 2 9) function corrected for 
changes in efficiency as a function of helicity angle. The 
BB backgrounds use an empirical non-parametric func- 
tion for AE, the masses and helicity angles. The contin- 
uum background tties shape is described by the function 
x-yl — x 2 exp[— £(1 — x 2 )] (with x = ntEs/Eg and £ a 
free parameter) [16| and a first- or third-order polyno- 
mial is used for AE and the helicity angles, respectively. 
The continuum invariant mass distributions contain real 
A"* candidates; we model the peaking mass component 
using the parameters extracted from the fit to the signal 
invariant mass distributions together with a second-order 
polynomial to represent the non-peaking component. 

We use the decay B° — ► D~ir + (D~ — ► K*°tt~) as a 
calibration channel to account for small differences be- 
tween MC simulation and reconstructed data. This de- 
cay has a similar topology to the modes under study and 
is selected using the same criteria as for K*°K*° but re- 
quiring the reconstructed A* invariant mass to be in 
the range 1.845 < mjf.o^ < 1.895 GeV/c 2 . We predict 
1860 ± 186 signal events and measure 1614 ± 47. 

We use MC-simulated events to study backgrounds 
from other B meson decays. The major charmless BB 
background to B° -> A* A* is B° -> (j>K*°, while 
charm BB backgrounds are effectively suppressed by the 
requirement that the two pions (and kaons) have oppo- 
site charge. For B° -> A* A* , B° -> (f>K*° remains the 
major charmless BB background, but a number of charm 
decays contaminate the signal, dominated by decays of 
the type B° -> D~K + and B~ -> D°K~. Given the 
uncertainty in the polarization and branching fractions 
of these backgrounds, we allow the BB background yield 
to float in the fit. 

A possible background is the decay B° — > 
A*°A*°(1430). We use the LASS parameterization for 
the A* (1430) lineshape, which consists of the A* (1430) 
resonance together with an effective-range non-resonant 
component [l7(. We apply the same selection criteria 
used for A* A* but require one of the A* candidates 
to have an invariant mass in the range 1.025 < ttik-k < 
1.53GeV/c 2 and perform an extended unbinned ML fit 
with the four variables m^s, AE, T , and the A* mass. 
We fit the LASS parameterization to the selected sig- 
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nal events in the i^*°(1430) mass range and extrapo- 
late to the K*° mass range. Interference effects between 
the K*° and the spin-0 final states (non-resonant and 
i^*°(1430)) integrate to zero as the acceptance of the de- 
tector and analysis is uniform. Assuming no interference, 
we expect 6 ± 5 B° -> K*°K*°(U30) events in the fit- 
ted B° — » K *°K*° signal region. The uncertainty on the 
contribution is calculated from the statistical error and 
the large uncertainty in the fitted LASS parameters used 
to describe the i^*°(1430) lincshape. We fix the yield in 
the final fit and vary the yield by its error to assess the 
systematic uncertainty. 

The continuum background PDF parameters that are 
allowed to vary are the T peak position, £ for toes, the 
slope of AE, and the polynomial coefficients and nor- 
malization describing the mass and hclicity angle distri- 
butions. We fit for B and /l directly and exploit the fact 
that B is less correlated with Jl than is cither the yield 
or efficiency taken separately. 

The total event sample consists of 7363 and 1390 events 
for B° -> K*°K*° and B° -> K*°K*°, respectively. The 
results of the ML fits are summarized in Tabled] The BB 
background yield agrees with the MC prediction within 
the statistical errors. The significance S of the signal 
is defined as S — 2Aln£, where Aln£ is the change 
in likelihood from the maximum value when the number 
of signal events is set to zero, corrected for the system- 
atic error defined below. The robustness of the signifi- 
cance estimate is cross-checked through fitting a series of 
toy MC ensembles generated from the fitted parameters. 
The significance of the B° — > x*°K*° branching fraction 
is 6<7, including statistical and systematic uncertainties. 
For B° -> K*°K*°, we compute the 90% C.L. upper 
limit as the branching fraction below which lies 90% of 
the total likelihood integral, taking into account the sys- 
tematic uncertainty. Figure [T] shows the projections of 
the fits onto toes, AE, K*° mass and cosine of the K*° 
helicity angle for B° -► K*°K*°. 

Systematic uncertainties in the branching fractions are 
dominated by our knowledge of the PDF modeling. Vary- 
ing the PDF parameters by their errors results in changes 
in the yields of 6.5% and 19.0% for B° -► K*°K*° and 
B° — ► K*°K*°, respectively. The largest contribution 
comes from the width of the K*°. 

The reconstruction efficiency depends on the decay po- 
larization. We calculate the efficiency using the measured 
polarization and assign a systematic error from the un- 
certainty on f L of 3.4% and 27.0% for B° -> K*°K*° 
and B° — > K*°K*°, respectively. Figure [2] shows the 
behavior of -2 In C(B, f L ) for the B° -> K*°K*° mode. 

The uncertainties in PDF modeling and /l are addi- 
tive in nature and affect the significance of the branch- 
ing fraction results. Multiplicative uncertainties in- 
clude reconstruction efficiency uncertainties from track- 
ing (3.2%) and particle identification (4.4%), track mul- 
tiplicity (1%), MC signal efficiency statistics (0.6%), and 
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K'° Mass (GeV/c 2 ) 

FIG. 1: (color online) Projections of the multidimensional 
fit onto (a) tties; (b) AE; (c) K*° mass; and (d) cosine of 
K*° helicity angle for B° — » K*° K*° events selected with a 
requirement on the signal-to-total likelihood probability ra- 
tio, optimized for each variable, with the plotted variable ex- 
cluded. The points with error bars show the data; the solid 
line shows signal-plus-background; the dashed line is the con- 
tinuum background; the hatched region is the signal; and the 
shaded region is the BB background. 
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FIG. 2: Distribution of -21n£(£,/ L ) for B° -> K*°K*° de- 
cays. The solid dot shows central value and the c urves give 
contours in steps of one sigma (A^/— 2 In C(B, /l) = 1). 



the number of BB pairs (1.1%). Variation of the ex- 
pected yield from B° — > K*° K*° (1430) events has a neg- 
ligible effect on the signal. 

The systematic uncertainty in fi is dominated by the 
PDF shape variations, which contribute 7% for B° — > 
K*°K*° and 20% for B Q -> K*°K*°. Other errors iden- 
tified above for the branching fraction have a very small 
effect on Jl and contribute in total 0.7%. The total sys- 
tematic error is summarized in Table HI 

In summary, we have measured the branching fraction 
B(B° -> K* Q K*°) = [1.28±g;|g(stat) ± O.ll(syst)] x lO" 6 
with a significance of 6a. We find the fraction of longi- 
tudinal polarization f L = 0.80±£^ (stat) ± 0.06 (syst). 
Both results are in agreement with the upper range of 
theoretical predictions. The 90% C.L. upper limit on the 



branching fraction B(B° ->■ K*°K*°) < 0.41 x 10~ 6 is 
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TABLE I: Summary of results: signal yield n s i g , the BB background yield n BB , signal reconstruction efficiency e (taking into 
account that B(K*° — > K + tt~) = 2/3), significance 5 (systematic uncertainties included), branching fraction £>, 90% C.L. 
upper limit for B° — ► K*°K*° branching fraction, and the longitudinal polarization Jl- The first error given is statistical and 
the second is systematic. 



Channel K K * K*"K* U 

~ g 33.5+g i 2.7 ±3.3 

n BB 1 9 ± 12 68 ± 29 

e (%) 6.8 6.4 

5 (a) 6 0.9 

B(10 _6 ) 1.28±g |g ± 0.11 O.lllg-l? ± 0.04 

UL 6(HT 6 ) - 0.41 

f L 0-80t°\l ± 0.06 1.0 ±1.0 



two orders of magnitude more stringent than previous 
measurements. 
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